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ABSTRACT 

We present a multi-frequency study of radio relics associated with the galaxy clusters A4038, A1664 
and A786. Radio images, integrated spectra, spectral index maps and fits to the integrated spectra in 
the framework of the adiabatic compression model are presented. Images of the relic in A4038 at 150, 
240 and 606 MHz with the GMRT have revealed extended ultra- steep spectrum (a ~ —1.8 to —2.7) 
emission of an extent 210 x 80 kpc 2 . The model of passively evolving radio lobes compressed by a 
shock best fits the integrated spectrum. The relic with a circular morphology at the outskirts of the 
cluster A1664 has an integrated spectral index of ~ — 1.10 ± 0.06 and is best fit by the model of radio 
lobes lurking for ~ 4 x 10 7 yr. The relic near A786 has a curved spectrum and is best fit by a model 
of radio lobes lurking for ~ 3 x 10 7 yr. At 4.7 GHz, a compact radio source, possibly the progenitor 
of the A786 relic, is detected near the center of the radio relic. The A786 radio relic is thus likely a 
lurking radio galaxy rather than a site of cosmological shock as has been considered in earlier studies. 
Subject headings: Galaxies: clusters: individual(A4038, A1664, A786)-Galaxies: clusters: intracluster 
medium- Radiation mechanisms: non-thermal- Radio continuum: general 



1. INTRODUCTION 

Diffuse extended radio sources with elongated or fil- 
amentary morphologies that cannot be identified with 
any unique optically detected galaxy are termed as ra- 
dio relics (e.g. see the review by Ferrari et al. 2008). 
The arclike radio relics that are typically found around 
massive and dynamically disturbed galaxy clusters (e.g., 
Bagchi et al. 2006; van Weeren et al. 2009) are proposed 
to be electrons accelerated at shock fronts by the diffu- 
sive shock acceleration mechanism (Ensslin et al. 1998). 
The accelerated electrons could be those from the ther- 
mal pool of the galaxy cluster or the old mildly rela- 
tivistic population due to active galaxies (Markevitch et 
al. 2005). Apart from the ~ Mpc size arc-like relics, 
there are relics that have smaller sizes (< 500 kpc) 
(e.g. Slee et al. 2001; van Weeren et al. 2009). These 
have been proposed to be lurking lobes of radio galax- 
ies or adiabatically compressed lobes of radio galaxies by 
shock waves (Ensslin & Gopal-Krishna 2001; Ensslin & 
Bruggen 2002). 

Regardless of their origin, radio relics are basically 
clouds of relativistic electrons and magnetic fields that 
evolve over timescales of ~ ten million years which 
is the lifetime of the synchrotron emitting electrons. 
Synchrotron emission, adiabatic expansion and inverse 
Compton scattering of the CMB photons are the mech- 
anisms by which the radio relics lose energy. The energy 
losses due to synchrotron and inverse-Compton scatter- 
ing depend on the square of the energy of the electrons 
and thus produce steepening in their spectra. There- 
fore to study the evolution of the radio relics, a multi- 
frequency approach is essential. In a multi-frequency 
study by Slee et al. (2001), integrated spectra of radio 
relics were modelled as evolving lobes of radio galax- 
ies losing energy by synchrotron emission and inverse- 
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Compton scattering (energy loss due to expansion was 
not considered). 

The spatial distribution of spectral indices over the ex- 
tent of the relic gives important clues about its origin. 
For example, a relic created at a shock front shows an 
edge of young, flat spectrum plasma and a gradual steep- 
ening of spectra behind the edge. Such trends have been 
mapped in a few arclike double radio relics such as A2345, 
A1240 (Bonafede et al. 2009), ZwCl 2341.1 + 0000, CIZA 
J2242.8+5301, ZwCl 0008.8 + 5215 (van Weeren et al. 
2009, 2010, 2011a) and a few arc-like single relics such 
as A2256 (Clarke & Ensslin 2006; Kale & Dwarakanath 
2010), A521 (Giacintucci et al. 2008) and A2744 (Orru 
et al. 2007). 

Spectral index distribution studies of single radio relics 
are rare (e.g., Cohen & Clarke 2011; Randall et al. 2010; 
van Weeren et al. 2011b). These relics show a wide va- 
riety in their morphologies and locations with respect 
to cluster centers. In this paper we present a multi- 
frequency study of three single radio relics which are 
associated with the galaxy clusters A4038, A1664 and 
A786. The relics and their host clusters are introduced 
in Sec. 1.1. The radio observations and data reduction 
are described in Sec. 2. Radio images of the relics are 
presented under Sec. 3. Integrated spectra are presented 
in Sec. 4 and model fits to the integrated spectra are 
described in Sec. 5. Results of each of the relic are dis- 
cussed in Sec. 6 and a summary is presented in Sec. 
7. We adopt a cosmology with Hq = 73 km s _1 Mpc _1 , 

^matter ~ 0.27 and ^vacuum ~ 0.73. 

1.1. The clusters and the relics 
The relic in A4038 (z= 0.03 (Struble and Rood 1999), 
K[,.i-2AkeV h £] ~ 2 -87 x 10 44 erg s" 1 (Reiprich & 
Bohringer 2002)) was discovered as an ultra-steep spec- 
trum source in surveys (160 MHz) by Slee & Siegman 
(1983). Based on further observations at 1.4 GHz (Slee 
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& Reynolds 1984; Slee & Roy 1998; Slee et al. 2001) 
(hereafter, SR98 and SL01), the extent of the relic was 
reported to be 56 kpc and a polarized intensity of 4% 
was detected. Spectral index (a, S oc v a ) over a narrow 
range of frequencies around 1.4 GHz was estimated to be 
-3.1. 

In A1664 (z= 0.1241 (Struble and Rood 1999), L x = 
5.36 x 10 44 erg s -1 (Pimbblet et al. 2006)) a relic was 
identified from the NVSS by Giovannini et al. (1999). It 
is known for its circular morphology (Govoni et al. 2001) 
unlike other relics which have elongated morphologies. It 
has not been imaged at any frequencies other than 1400 
MHz and thus the spectral index is unknown. 

The relic near A786 was first identified and studied by 
Dewdney et al. (1991) and Harris et al. (1993). This relic 
is located at a distance of ~ 5 Mpc from the center of 
the cluster A786 (z= 0.128 (Struble and Rood 1999), log 
(L x ) ~ 44.8). The relic has been discussed as a site of a 
large scale cosmological shock (Ensslin et al. 1998). This 
relic is also referred to as '0917+75' in literature. 

2. RADIO OBSERVATIONS AND DATA REDUCTION 

The Giant Meterwave Radio Telescope (GMRT) was 
used to carry out observations of the relics in the fre- 
quency bands of 150, 240, 325, 606 and 1288 MHz. The 
shortest baseline of 100 m and the longest baseline of 
25 km of the GMRT in a hybrid configuration provide 
an opportunity to image extended sources such as radio 
relics with good resolution and sensitivity at multiple 
frequencies. The Westerbork Synthesis Radio Telescope 
(WSRT) was used to carry out observations of the relic 
in A786 at 345 MHz. Apart from these, archival data 
from the Very Large Array (VLA) at 4.7 GHz and im- 
ages from the NRAO VLA Sky Survey (1.4 GHz, NVSS, 
Condon et al. (1998)) and the VLA Low Frequency Sky 
Survey (74 MHz, VLSS, Cohen et al. (2007)) were used. 
A summary is given in Table 1. 

Reduction of all the data were carried out using NRAO 
AIPS (Astronomical Image Processing System). Stan- 
dard steps of data editing and complex gain calibra- 
tion were carried out. The data lost due to radio fre- 
quency interference (RFI) varied between 5 — 35% over 
the range of frequencies 1288 to 150 MHz at the GMRT. 
The WSRT observation of the relic in A786 consisted of 
8 bands, each of band width 10 MHz, between 310 and 
390 MHz. One frequency band of 10 MHz had to be 
discarded due to RFI. The latest VLA values were used 
for setting flux calibration (AIPS version 31DEC10). In 
the frequency bands, 150, 240, 325 and 606 MHz, wide 
field imaging was carried out to overcome the effects of 
non-coplanarity. The visibilities were weighted accord- 
ing to the 'ROBUST=0' scheme implemented in AIPS. 
Self-calibration was carried out to obtain deeper images 
at all the frequencies. 

The angular extents of the relics (~ 5', 6' and 10', see 
Sec. 3) are much smaller as compared to the largest 
angular extents that can be fully imaged with the GMRT 
at 240 (~ 44'), 325(~ 32') and 606 (~ 17') MHz, with 
the WSRT at 345 (~ 95') MHz and with the VLA-D at 
1400 (~ 15') MHz. Therefore no flux is missed due to 
missing short baselines. Images at multiple frequencies 
made with the same resolution, corrected for the primary 
beam gains and blanked in regions with flux densities 
below 5a were used to construct spectral index maps. 



3. RADIO IMAGES AND SPECTRAL INDEX MAPS 

3.1. Abell 4038 

Images of A4038 using the GMRT at 150 and 1288 
MHz are presented in Fig. 1 as contours and grey-scale 
respectively. The elongated, filamentary, 'hook' shaped 
structure detected at 1288 MHz is the radio relic in 
A4038 as reported earlier in the literature (SR98). At 
150 MHz extended emission about thrice the size of that 
at 1288 MHz is detected. The extent of the relic de- 
tected at 150 MHz is 133 kpc (~ 3.8') along N-S and 
80.5 kpc (~ 2.3') along E-W. The positions of discrete 
radio sources (unresolved sources or those identified with 
optically detected galaxies) are marked by crosses. The 
northern cross marks a compact radio source which does 
not have any optical counterpart (A4038_10, SL01) and is 
presumed to be a background source unrelated to A4038. 
The southern cross marks a radio source identified with 
the cD galaxy in A4038 (A4038_ll, SL01). In Fig. 2 
(left), shown in grey-scale is the DSS R band image of 
the A4038 field. The cD galaxy can be clearly seen in 
the R band image. 

In the high resolution image at 240 MHz (Fig. 2, left) 
blobs of emission are seen scattered to the NW of the 
cD (A4038_ll) which indicate the presence of extended 
features. We produced low resolution (25" x 25") im- 
ages at 240 MHz (Fig. 2, right) to confirm the presence 
of the suspected extended features. We refer to the low 
surface brightness extension toward north west as the 
'NW-Extension' or 'NW-Ext' for brevity, hereafter. The 
NW-Ext has an extent of 84x56 kpc 2 . The average sur- 
face brightness of the NW-Ext at 240 MHz is -7 mJy 
beam" 1 (beam= 25" x 25"). The image at 606 MHz 
convolved to a resolution of 25" x 25" is presented in 
Fig. 3 (contours). The high resolution 606 MHz image 
is shown in grey-scale. The low surface brightness NW- 
Ext is not detected at 606 and 150 MHz. The surface 
brightness detections limited to ~ 2 mJy beam -1 at 606 
MHz and to ~ 20 mJy beam -1 at 150 MHz resulted in 
non-detections of the NW-Ext at these frequencies. The 
non-detection of NW-Ext at 606 MHz implies an upper- 
limit on the spectral index of NW-Ext of a^o < — 1-3. 
The NVSS (Condon et al. 1998) and the VLSS (Cohen 
et al. 2007) images at 1400 and 74 MHz respectively, 
are presented in Fig. 4. Extended emission is clearly 
detected in both the images. 

We refer to the high surface brightness regions of the 
relic detected at the frequencies 150, 240 and 606 MHz 
(Figs. 1, 2 and 3) as the 'main relic'. A reference to 'the 
relic', hereafter will mean the 'main relic' together with 
the NW-Ext. 

The spectral index map of A4038 between 606 and 240 
MHz with a resolution of 25" x 25" and the corresponding 
uncertainty map are presented in Fig. 5. The positions 
of the discrete sources are marked by '+'s and are not 
subtracted from the images. The map shows a complex 
distribution of spectral indices. The central region of the 
relic has a spectral index of -1.3 and it steepens toward 
north and south reaching spectral indices of -2.7. 

3.2. Abell 1664 

Images of A1664 at 150 and 325 MHz using the GMRT 
and at 1400 MHz from the NVSS are presented in Fig. 6. 
The relic has a circular morphology and has an angular 
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extent of ~ 6' (linear extent ~ 0.927 Mpc). To the north 
east, the relic appears connected to an unresolved radio 
source associated with the galaxy A1664.3 (Pimbblet et 
al. 2006). The high resolution image at 325 MHz with 
the GMRT resolves this compact source (Fig. 7). The 
spectral index map between 325 and 1400 MHz and the 
corresponding uncertainty map are given in Fig. 8. The 
relic has a flat spectral index distribution (~ —0.5 to 
— 1.0) with small pockets of steep spectral index. The 
source A1664_3 has a spectral index of -0.5. 

3.3. Abell 786 

The images of the relic associated with A786 at 150, 
345, 606 and 1400 MHz are presented in Fig. 9. The 
largest extent of the relic is detected at 345 MHz (~ 10'). 
A central region with higher surface brightness is seen at 
345 and 1400 MHz with low surface brightness extensions 
along east and west. A large fraction (~ 40%) of short 
baseline data at 150 MHz was corrupted due to RFI and 
had to be discarded. The sensitivity in the 150 MHz 
image is ~ 10 mJy beam -1 . The bright central portion 
of the relic is detected but the weak extensions along 
the east and west are not detected at 150 MHz. To the 
southwest of the relic is a bright unresolved radio source 
(source R5, Harris et al. (1993)). It is presumably a 
background source not related to the relic. 

A high resolution image at 606 MHz (7" x 4", PA 
= 0.16°) is shown in Fig. 10 (left). Two unresolved 
sources (marked as 'twin sources' in Fig. 10 and by 
in Fig. 9) were detected at the northern edge 
of the relic. The positions (J2000) of these sources 
are RA 09h22ml9.096s DEC ±75d01'41.95" and RA 
09h22m24.77s DEC ±75d01'43.92" (separation - 20"). 
No counterparts to these are found in the DSS and 
2MASS images. We also produced an image at 4.7 GHz 
using archival data from the VLA (D-array) (Fig. 10, 
right). A compact source is detected in the central re- 
gion of the relic (RA 09h22m09.27s DEC 74d59 / 48.99", 
J2000) with extended emission of linear extent ~ 350 
kpc surrounding it. A counterpart to this source is de- 
tected in 2MASS (total magnitudes in J, H and K s bands 
are 13.976, 13.584 and 12.701, respectively), however no 
information on the redshift of it is available. From the 
K s — z correlation (de Breuck et al. 2002), the K s band 
magnitude implies a redshift of ~ 0.1 (redshift of A786 is 
0.128). Spectral index maps between 345 and 606 MHz 
and between 606 and 1400 MHz were produced and are 
shown in Figs. 11 and 12 respectively. The distribution 
of spectral indices over the relic is complex. 

4. INTEGRATED SPECTRA 
4.1. A4038 

The total flux density in a box enclosing the relic and 
two discrete sources (A4038_ll and A4038_10) in A4038 
was found using the images with resolutions 25" x 25" at 
each of the frequencies (150, 240 and 606 MHz). Total 
flux densities were also estimated at 1400 and 74 MHz 
using the NVSS and the VLSS images, respectively, using 
boxes of same dimensions. To obtain the flux density of 
only the relic, contribution from the sources A4038_ll 
and A4038_10 had to be subtracted from the total flux 
density in the box. The high resolution images at 1288 
and 606 MHz were used to find the flux densities of the 



sources A4038_ll and A4038_10. Since these sources are 
embedded in the diffuse emission of the relic at 74, 150, 
240 and 1400 MHz and could not be resolved, the spectral 
index between 1288 and 606 MHz was used to estimate 
their contributions at the other frequencies. The spectral 
indices (a^of) of A4038_ll and A4038_10 were found to 
be -0.55 and -0.88, respectively. It was assumed that the 
spectral indices are constant over the frequency range 
of 74 to 1400 MHz. The flux densities at 29, 327, 408 
and 843 MHz corrected for the contributions from the 
sources A4038_ll and A4038_10 and also other sources 
in the beams as reported by SL01 are given in Table 2. 
The integrated spectrum of the relic is plotted in Fig. 13 
(points). The power of the relic at 1.4 GHz is 1.13 x 10 23 
W Hz" 1 . 

The flux densities of the relic at 80 and 160 MHz are 
reported to be 19.0 ±2.7 Jy and 4.3 ±0.5 Jy respectively, 
by SL01. However since we have better estimates at 74 
and 150 MHz (Table 2), we do not use the SL01 values 
at 80 and 160 MHz. 

4.2. A1664 

The flux density in a box enclosing the relic and the 
radio source A1664_3 was estimated. From higher reso- 
lution images, the flux density of A1664_3 was estimated 
and subtracted. These flux densities are listed in Table 
3. The integrated spectrum of the relic is plotted in Fig. 
15 (points). The spectral indices of the relic, af 2 Q and 

^325° are - 1 - 32 ± °- 03 and -°- 98 ± °- 01 > respectively. 
The power of the relic at 1.4 GHz is 3.41 x 10 24 W Hz" 1 . 
In the highest resolution map (12" x 6", rms 0.6 mJy 
beam -1 ) at 325 MHz, an unresolved source is detected 
in the region of the relic (marked in Fig. 6). The po- 
sition (J2000) of the source is RA 13h03m30.379s DEC 
— 24d21 / 54.29" and the peak and the integrated flux den- 
sities are 6.14 ± 0.06 and 8.1 ± 1.3 mJy beam -1 , respec- 
tively (about 2% of the total flux density of the relic). 
This source does not have a detectable counterpart in 
the DSS R band image and in the 2MASS J, H and K s 
band images. 

4.3. A786 

The flux density in a box enclosing the relic, the twin 
discrete sources and the central source detected at 4.7 
GHz was determined from the images at 150, 345, 606 
and 1400 MHz. The total flux densities of the twin dis- 
crete sources in the relic region detected at 606 MHz and 
4.7 GHz are 7.6 mJy and 1.5 mJy, respectively. Using the 
spectral index of -0.78 between 606 MHz and 4.7 GHz of 
the twin sources, the flux densities at 150, 345 and 1400 
MHz were estimated and subtracted from the flux densi- 
ties of the relic. The integrated flux densities of the relic 
after this correction are reported in Table 3. The power 
of the relic at 1.4 GHz is 3.59 x 10 24 W Hz -1 . The source 
detected at 4.7 GHz in the central region of the relic has 
a peak flux density of 0.34 ±0.06 mJy beam -1 and a total 
(including the surrounding extended emission) flux den- 
sity of 7.2 ±0.3 mJy beam -1 . Since the full extent of the 
relic is not detected at 4.7 GHz we do not report total 
flux density of the relic at 4.7 GHz. The central compact 
source is not detected at lower frequencies, therefore no 
correction due to this source is made to the flux densities 
of the relic reported at lower frequencies. The integrated 
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spectrum of the relic is plotted in Fig. 16 (points). 
5. ADIABATIC COMPRESSION MODEL FITS 

Radio relics are proposed to be either remnants of radio 
galaxies or sites of accretion/ merger shocks where parti- 
cles are accelerated to relativist ic energies. The remnants 
of radio galaxies can be ageing lobes or lobes compressed 
by shocks in the intracluster medium (ICM). In the case 
of A4038, SL01 attempted a fit considering the relic as a 
remnant of a radio galaxy. The SL01 fits were based on 
estimates of volume and magnetic field using the extent 
of the relic (LLS=56 kpc) as detected by SL01 at 1.4 
GHz. This extent of the relic is an underestimate due to 
the discoveries of extended features of the relic at lower 
frequencies. Further the models considered by SL01 ne- 
glected energy losses due to expansion of the lobes. We 
describe our attempt to fit the integrated spectrum of the 
relic with a model which overcomes the above mentioned 
drawbacks. 

Revival of fossil radio lobes by adiabatic compression 
by passage of shocks is one of the models favored for radio 
relics (Ensslin & Gopal-Krishna 2001) (EG01 hereafter). 
The framework of EG01 considers revival of fossil radio 
cocoon by adiabatic compression due to the passage of 
shocks in the ICM. The evolution of the spectrum of a ra- 
dio galaxy before and after the jets cease to be active and 
after adiabatic compression of the lobes can be obtained. 
Energy losses due to synchrotron, inverse Compton (IC) 
and expansion are accounted for. There are five phases 
of evolution, viz. injection, expansion, lurking, flashing 
and fading. The salient features of the model and the 
application in the context of radio relics can be found in 
Kale & Dwarakanath (2009). 

The largest linear extent of the relic in A4038 is 210 
kpc. It is is located close to the cD galaxy and resides 
almost at the cluster center (in projection). The volume 
of the complete relic was approximated to be the sum 
of the volumes of the main relic (110 x 75 x 75 kpc 3 ) 
and that of the NW-Ext (25 3 kpc 3 ). A magnetic field of 
3/iG was estimated using the minimum energy condition 
under standard assumptions. 

Best fits to the integrated spectrum of A4038 relic in 
'flashing' and 'fading' phases were obtained. The reduced 
Chi -square (Xred) was estimated to be 5.96 and 5.49 in 
the two fits respectively. The parameters used in the 
model for the best fit are listed in Table 4. The model 
fits along with the observed spectrum are plotted in Fig. 
13 and 14 (left). The best fits are denoted by black lines. 
The model spectra in other phases are also plotted in the 
same figures for comparison. The right hand side panels 
in Figs. 13 and 14 show the best fits along with model 
spectra obtained when the values of the duration of the 
phases were changed by ±30%. Note the deviation from 
the observed spectrum. 

The relic in A1664 is at a distance of ~ 1.13 Mpc from 
the center of the cluster. Best fit to the spectrum was 
found when the relic was assumed to be in the 'lurking' 
phase (Fig. 15). The best fit parameters are listed in 
Table 4. The duration of the 'lurking' phase was found 
to be 4 x 10 7 yr. 

The EG01 model was used to fit the integrated spec- 
trum of the relic near A786. A best fit was obtained 
when the relic was assumed to be in the 'lurking' phase. 
The best fit parameters are reported in Table 4 and the 



fit is shown in Fig. 16 (left). The plots showing devi- 
ation from the best fit when the timescale was changed 
by ±30% are shown in Fig. 16 (right). The duration of 
the lurking phase phase was found to be 6 x 10 7 yr. 

6. DISCUSSION 

The imaging of the three radio relics in galaxy clusters 
A4038, A1664 and A786 at low frequencies (< 1 GHz) 
has revealed several new properties of the relics. Im- 
ages of comparable sensitivity and resolution were pro- 
duced by the use of the GMRT, the WSRT and the 
VLA over the frequency range of 150-1400 MHz. The 
integrated spectra sampled at more than three frequen- 
cies were produced. The spectra are curved and steep 
(a < —1) and unlike a flat power-law expected of young 
freshly accelerated plasma (a ~ —0.7). The adiabatic 
compression model was used to fit the integrated spec- 
tra. The models of radio galaxy cocoons of ages > 10 7 
yr fit the spectra. The spectral index maps show com- 
plex distributions. The gradients in the spectral index 
maps do not show edges of flat spectral indices that are 
characteristic of plasma accelerated at shock fronts. The 
spectral indices are steep and vary across the extent of 
the relic by as much as 1 (Aa ~ 1). These properties im- 
ply that these relics contain aged relativistic plasma and 
can be termed as 'relics of radio galaxies'. The model 
of shock accelerated plasma is disfavored by the observa- 
tions. The properties of the relics are discussed below. 

6.1. Abell 4038 

Our sensitive low frequency observations with the 
GMRT have led to the discovery of steep spectrum emis- 
sion in addition to the relic known from earlier high fre- 
quency studies. The largest linear extent of the relic in 
A4038 as seen at 240 MHz is ~ 210 kpc as opposed to 56 
kpc detected at 1.4 GHz. The relic is elongated along the 
N-S direction. A low surface brightness extension to the 
relic toward the NW has been discovered at 240 MHz. 
Due to limitations of the sensitivities of observations at 
150 and 606 MHz, the NW-Ext was not detected at these 
frequencies. Deeper observations at these frequencies are 
required to detect and determine the spectral index of the 
NW-Ext. 

We used the adiabatic compression model (EG01) to 
fit the integrated spectrum of the relic. The EG01 model 
accounts for the energy losses due to expansion and also 
the resulting change in the magnetic field due to change 
in volume. Further, it also considers the case in which 
a shock in the ICM compresses lurking radio cocoons 
and leads to enhanced emission from them. The best fits 
were obtained in the 'flashing' and the 'fading' phases. 
The durations of the phases of injection, expansion and 
lurking are the same in both the cases. However the 
extent of the compression in the 'flashing' phases dif- 
fers. In the case of 'flashing' phase being the best fit, 
the volume changes by a factor of ~ 9.3, whereas in the 
case of 'fading' being the best fit, the change is by a 
factor of ~ 3.5. The relation between the compression 
factor (C) and the pressure ratio (P1/P2) is given by, 
C = (Pi/i^) 3 ^ 4 - This implies that the pressure jumps 
at the shock that caused the compression are factors of 
~ 20 and ~ 5 in the 'flashing' and the 'fading' best fit 
cases, respectively. Using the pressure ratios and the fac- 
tor 7 = 5/3 (adiabatic index for thermal gas), the Mach 
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number M (Eqn. 22, EG01) was found to be ~ 4 and ~ 2 
in the cases of 'flashing' and 'fading' best fit respectively. 
It should be noted that the role of shock in this context is 
to compress the cocoon. The speed of sound in the radio 
cocoon being very high(~ c/y/3), the relativistic plasma 
is not shocked, but is compressed as a result of an en- 
counter with a shock (EG01). No direct acceleration at 
the shock front is believed to be operative. The change 
in the spectral luminosity by factors of 10 is achieved in 
the compression phase making the relic detectable over 
the observed frequency range. 

The remnants of radio galaxies in clusters have been 
found in the form of bubbles in the ICM (eg. McNa- 
mara (2002)). The jets of radio galaxies inflate the bub- 
bles in the ICM which are detected in high resolution 
X-ray images as cavities. The cavities (bubbles) are of- 
ten found to be filled with relativistic plasma detectable 
at radio wavelengths. These bubbles can get detached 
from the place where they formed and rise in the ICM 
due to buoyancy. The main relic in A4038 could be such 
a bubble and the NW-Ext could be an older bubble that 
has risen by buoyancy. In a study of radio bubbles in 
galaxy clusters, A4038 has been classified as a source 
with a central radio source but no detected cavities in 
X-ray observations with Chandra (Dunn et al. 2005). In 
a later study it has been remarked that the sensitivity 
of the X-ray observations may have been insufficient in 
revealing the cavities since no extended radio emission 
was present to provide clues about the location of cavi- 
ties (Dunn & Fabian 2006). Given that we have detected 
extended steep spectrum radio emission in A4038, a re- 
analysis of the existing X-ray observations or a deeper 
observation may lead to detection of cavities. 

The relativistic plasma that forms the relics has been 
considered to originate from a radio galaxy in the EG01 
model. Over time, the host galaxy would have moved 
and may not have an AGN any more, preventing iden- 
tification. The radio images of A4038 relic indicate a 
diffuse bridge connecting the relic to the cD. The bridge 
cannot be an active jet considering that it has a steep 
spectral index and is not detectable at higher frequen- 
cies (young relativistic plasma in the jet would have a 
flat spectrum). The cD galaxy could be the progenitor 
of the relic but the evidence is not sufficient to rule out 
other progenitors. 

6.2. Abell 1664 

The 150 and 325 MHz images presented above are the 
first low frequency images of the relic in A1664. These 
have yielded an estimate of the spectral index of the 
relic between 150 and 1400 MHz - -1.10 ± 0.06. We 
attempted a model fit to the spectrum using the EG01 
model. The best fit found in the 'lurking' phase implies 
that no shock compression is required to produce the 
relic. This is consistent with the fact that there has not 
been any recent dynamical activity in the cluster as also 
implied by the cooling flow at the center of A1664 (Allen 
et al. 1995). The cluster A85 which is also a cooling flow 
cluster hosts a peripheral relic. The main difference be- 
tween A85 and A1664 relics is that the relic in A85 has 
an extent of ~ 100 kpc whereas A1664 relic has an extent 
> 0.92 Mpc. 

The NVSS image overlaid on ROSAT (0.1-2.4 keV) im- 
age of A1664 is shown in Fig. 17. The cluster center is 



to the north of the relic. The X-ray surface brightness in 
the cluster center (Chandra image) is elongated and can 
be approximated to be an ellipse, the major axis of which 
is at a P. A. of 25° (Kirkpatrick et al. 2009). The relic 
is situated along the same axis toward SW at a distance 
of 1.13 Mpc from the cluster center. The elongation in 
X-rays is consistent with the optical major axis of the 
galaxy at the cluster center and thus the elongation may 
be due to processes local to the cD or a non relaxed state 
(Pierre & Starck 1998). The location of the relic at the 
outskirts of the cluster and close to a radio source is sim- 
ilar to the relic in Coma cluster. The Coma cluster relic 
is at a distance of 2.2 Mpc from the cluster center and 
is possibly associated with a radio galaxy seen close to 
the relic (Ensslin et al. 1998; EG01). The cosmologi- 
cal structure formation shocks are believed to be found 
at cluster peripheries. These could compress cocoons of 
radio galaxies leading to the formation of relics and/or 
could accelerate the particles in the thermal gas to form 
relics. The Coma relic is elongated and has been pro- 
posed to have been compressed by an accretion shock 
(Ensslin et al. 1998). The A1664 relic shows similarity 
in terms of its location relative to the cluster center to 
the Coma relic, but it is morphologically different and 
may have a completely different origin. 

The relic appears connected to the galaxy A1664_3 
in the images. However, no evidence for jets from the 
galaxy connecting to the relic has been found. It is pos- 
sible that the A1664 relic is a blob of radio emission from 
the past activity of an AGN in the cluster. Deeper X-ray 
studies in the region of the relic are required to detect if 
there is a cavity corresponding to the A1664 relic. The 
model fit to the spectrum in the EG01 framework and 
the morphology of the A1664 imply that the relic is a 
lurking lobe of a radio galaxy. 

6.3. Abell 786 

The relic is about 5 Mpc from the center of the cluster 
A786 in the plane of the sky. In Fig. 18 the 345 MHz 
contours are overlaid on the ROSAT HRI image (0.1-2.4 
keV). There are two optical galaxies near the brightest 
central region of the relic, called D and E by Harris et al. 
(1993), which have redshifts consistent with that of the 
cluster A786. For this reason the relic has been consid- 
ered associated with the cluster A786. 

The integrated spectrum of the relic near A786 is 
curved. The spectral index between 150 and 606 MHz is 
~ —0.80 and it steepens to a value of ~ —1.27 between 
606 and 1400 MHz. A best fit in 'lurking' phase was 
obtained. According to the model fit the radio galaxy 
responsible for depositing the relativistic plasma in the 
ICM had been active for ^ 1.5 x 10 7 yr. The total time 
elapsed since the jets of the radio galaxy ceased to be ac- 
tive is ~ 5 x 10 T yr. Within this timescale the progenitor 
galaxy would have travelled a distance of 50 kpc with a 
velocity of 1000 km s _1 (typical at cluster center) and 5 
kpc with a velocity of 100 km s _1 (typical velocities in 
galaxy groups). Since at a distance of 5 Mpc from the 
cluster center the velocities of galaxies are low (similar to 
galaxy groups), it is likely that the progenitor lies within 
the extent of the relic. The compact radio source with 
extended features surrounding it, detected at 4.7 GHz 
could be the progenitor AGN and the relic, a fossil lobe 
associated with it. 
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The relic near A786 is one of the largest diffuse emis- 
sion (LLS~ 1.6 Mpc) regions known in the Universe. 
Such large scale filamentary radio emission is rare and 
has been associated with filaments of galaxies or cluster 
peripheries (eg. ZwCl 2341.1+0000, Bagchi et al. 2002 
and A3376, Bagchi et al. 2006). Structure formation 
shocks are believed to be responsible for such sources. 
Ensslin et al. (1998) have discussed this relic as a candi- 
date for accretion shock compressed relic. However due 
to lack of information regarding the temperature and ve- 
locity dispersion of the cluster A786, the expected shock 
radius and polarization could not be determined (Ensslin 
et al. 1998). The spectral index maps of the relic pre- 
sented here do not show gradients as would be expected 
in the case of shock acceleration. Also the integrated 
spectrum of the relic is consistent with that of a 'lurk- 
ing' radio galaxy lobe. Thus, it seems that this relic is 
not a site of cosmological structure formation shock ac- 
celeration. 

7. SUMMARY 

A multi-frequency study of the three radio relics asso- 
ciated with the galaxy clusters, A4038, A1664 and A786 
is presented. Radio images and spectral index maps con- 
structed from those are presented and discussed. We 
have used the adiabatic compression model to fit the in- 
tegrated spectra of the relics. 

The 150, 240 and 606 MHz images of A4038 have revo- 
lutionized the picture of the relic known earlier. New fea- 
tures of the relic are uncovered extending the largest de- 
tected linear size of the relic from 56 kpc to 210 kpc. The 
newly discovered emission has spectral indices a < —1.8. 
The model of radio lobes compressed by shock best fits 
the integrated spectrum of A4038 relic. The relic in 
A1664 was imaged at 150 and 325 MHz with the GMRT. 
These images and the NVSS image at 1400 MHz were 
used to construct the integrated spectrum of the relic in 
A1664. A best fit to the integrated spectrum of the relic 
in A1664 was obtained in the 'lurking' phase. The im- 
ages of the relic near A786 at 150, 345, 606 and 1400 MHz 
were used to construct the integrated spectrum. The in- 
tegrated spectrum was best fit in the EG01 framework 
in the 'lurking' phase. A source detected at 4.7 GHz in 
the central region of the relic that has an optically iden- 
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TABLE 2 

Integrated flux densities of the relic in A4038. 



Freq. (MHz) 
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TABLE 3 

Integrated flux densities of the relics in A1664 and in A786. 
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TABLE 4 

Adiabatic compression model fit parameters for the relics. 
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Fig. 1. A4038: GMRT 150 MHz (25" x 25") contours overlaid on GMRT 1288 MHz (3.8" x 2.2", Position angle= 36°) image shown 
grey- scale. The contours are at -15, 15, 20, 25, 40, 60, 120, 240, 480, 770 mJy beam -1 . Grey-scale extends from 0.25 to 0.8 mJy beam - 
The plus signs mark the positions of discrete radio sources A4038_10 and A4038_ll. 
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Fig. 2. A4038: Left: GMRT 240 MHz (15" x 10", Position angle= -11°) contours overlaid on optical R band image (DSS) shown in 
grey- scale. The contours are at -2.7, 2.7, 3.6, 4.5, 7.2, 10.8, 21.6, 43.2, 86.4, 138.6 mJy beam- 1 . The grey-scale extends from 2500 to 
10000 counts. Right: GMRT 240 MHz image (25" x 25"). The contours are at -4.5, 4.5, 6.0, 7.5, 12.0, 18.0, 36.0,72.0, 144.0, 231.0 mJy 
beam -1 . 
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Fig. 3. A4038: GMRT 606 MHz (25" x 25") low resolution image shown in contours overlaid on GMRT 606 MHz (6.6" x 3.5", Position 
angle= 28°) high resolution image in grey-scale. The contours are at -1.2, 1.2, 1.6, 2.0, 3.2, 4.8, 9.6, 19.2, 38.4, 61.6 mJy beam -1 . The 
grey-scale extends over the range of 0.75 to 3 mJy beam -1 . 
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Fig. 4. A4038: Left: Image of A4038 from the NVSS (1400 MHz). The synthesized beam is 45" x 45". Right: VLSS (74 MHz) image 
of A4038. The synthesized beam is 80" x 80". Contours are at (left) -1.35, 1.35, 1.80, 2.25, 3.6, 5.4, 10.8, 21.6 mJy beam -1 and at (right) 
-300, 300, 400, 500, 800, 1200, 2400, 4800 mJy beam" 1 . 
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Fig. 5. A4038: Left: Spectral index image of A4038 between 606 and 240 MHz. The synthesized beam is 25" x 25". Right: Spectral 
index uncertainty image. 
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RIGHT ASCENSION (J2000) RIGHT ASCENSION (J2000) RIGHT ASCENSION (J2000) 

Fig. 6. A1664: Left: GMRT 150 MHz contours. Contour levels are at -21, 21, 42, 63, 84, 105, 210, 420, 840 mJy beam -1 (synthesized 
beam = 54" x 54"). Center: GMRT 325 MHz contour levels at -4.5, 4.5, 9.0, 13.5, 18.0, 22.5, 45.0, 90.0, 180.0 mJy beam" 1 (synthesized 
beam = 54" x 54"). Grey-scale extends from 3 to 40 mJy beam -1 . The location of the unresolved source detected in the high resolution 
image at 325 MHz is marked ('+'). Right: NVSS 1400 MHz contour levels at -1.05, 1.05, 2.10, 3.15, 4.20, 5.25, 10.50, 21.00, 42.00 mJy 
beam - 1 (synthesized beam = 45" x 45"). The source A1664_3 is seen at the north of the relic. 
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Fig. 7. — A1664: The source A1664_3 as seen at 325 MHz with high resolution (synthesized beam= 12" x 6", Position angle= 5°) . The 
325 MHz contours (-3, 3, 6, 9, 15, 20, 25, 30, 30, 35, 40 mJy beam -1 ) are overlaid on the DSS R band image of the same field shown in 
grey-scale. 
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Fig. 8. A1664: Left: Spectral index image of A1664 between 1400 and 325 MHz. The synthesized beam is 25" x 25". Right: Spectral 
index uncertainty image. 
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Fig. 9. A786: Top Left: GMRT 150 MHz. Contour levels are at -30, 30, 60, 90, 120, 150, 200, 400, 600, 800, 1600 mJy 
beam- 1 (synthesized beam= 67" x 52", P.A.= 22°). Top Right: WSRT 345 MHz. Contour levels are at -2.4, 2.4, 4.8, 7.2, 9.6, 12.0, 
16.0, 32.0, 48.0, 64.0, 128.0 mJy beam -1 (synthesized beam= 64" x 64"). Bottom Left: GMRT 606 MHz contours overlaid on DSS R band 
image in grey-scale. Contour levels are at -6.9, 6.9, 20.7, 27.6, 34.5, 46.0, 92.0, 184.0, 368.0 mJy beam" 1 (synthesized beam= 64" x 64"). 
Bottom Right: VLA 1400 MHz. Contour levels are at -0.3, 0.3, 0.6, 0.9, 1.2, 1.8, 2.4, 3.0, 4.0, 5.0, 7.0, 10.0, 12.0, 20.0, 50.0 mJy 
beam -1 (synthesized beam = 64" x 64"). 
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Fig. 10. — A786: Images of the region of the relic in A786 at 606 MHz (left) and 4.7 GHz (right). A central radio source with extended 
emission surrounding it is detected at 4.7 GHz. The twin sources are detected at 606 MHz but the central source is not. The synthesized 
beams are 7" x 4" (P. A.= 0.16°) and 11" x 11" (P. A.= -64°) in the left and in the right panels respectively. 
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Fig. 12. A786: Left: Spectral index image of A786 between 606 and 1400 MHz. The synthesized beam is 64" x 64". Right: Spectral 
index uncertainty image. 
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Fig. 13. — A4038: Left: The best fit model spectrum in 'flashing' phase along with spectra in other phases. Right: The observed 
integrated spectrum of the relic with the best fit and the model spectra with A^3 db 30%. The red points mark values estimated using our 
GMRT images (150, 240 and 606 MHz) and the images from the VLSS (74 MHz) and the NVSS (1400 MHz). Black points mark values 
from Slee et al (2001). See the electronic edition of the Journal for a color version of this figure. 




Fig. 14. — A4038: Left: The best fit model spectrum in 'fading' phase along with spectra in other phases. Right: The observed integrated 
spectrum of the relic with the best fit and the model spectra with At 4 ± 30%. See the electronic edition of the Journal for a color version 
of this figure. 
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Fig. 15. — A1664: Left: The integrated spectrum of the relic in A1664 is plotted (Table 3). The lines are the spectra obtained from the 
EG01 model. Solid black line is the best fit ('lurking Phase') to the integrated spectrum of the relic in A664. Other lines are model spectra 
in the other phases of evolution. Right: The best fit model spectrum (black line) with the model spectra obtained when the duration of 
the best fit phase (A^) was varied by ±30%. The integrated spectrum (points) is also plotted for reference. See the electronic edition of 
the Journal for a color version of this figure. 
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Fig. 16. — A786: Integrated spectrum of the relic near A786 is shown by points with error bars. Left Each curve is a spectrum in a Phase 
(labelled in the legend) in the model. The solid black line is the best fit spectrum in the 'lurking' phase. Right The curves representing 
spectra when the duration of the best fit phase was changed by ±30%. The best fit curve and the integrated spectrum are plotted for 
reference. See the electronic edition of the Journal for a color version of this figure. 
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Fig. 17. — A1664: 1400 MHz image (contours) overlaid on ROSAT (0.1-2.4 keV) image in grey-scale. Contour levels are at -1.05, 1.05, 
2.10, 3.15, 4.20, 5.25, 10.50, 21.00, 42.00 mJy beam" 1 . 
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Fig. 18. A786: Contours at 345 MHz (WSRT) overlaid on the ROSAT HRI (0.1-2.4 keV) image. The cluster A786 lies to the southeast 
of the relic and is seen as a bright source in X-rays. A radio galaxy resides at the center of A786. Contour levels are at -2.4, 2.4, 4.8, 7.2, 
9.6, 12.0, 14.4, 19.2, 24.0, 32.0, 40.0, 64.0, 128.0, 256.0, 512.0 mJy beam" 1 (synthesized beam= 52" x 51", P.A.= 7°) 



